Anisotropies in the Cosmic Neutrino Background after WMAP 5-year Data 
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We search for the presence of cosmological neutrino background (CNB) anisotropies in recent 
WMAP 5-year data using their signature imprinted on modifications to cosmic microwave back- 
ground (CMB) anisotropy power spectrum. By parametrizing the neutrino background anisotropies 
with the speed viscosity parameter c v i s , we find that the WMAP 5-year data alone provide only 
a weak indication for CNB anisotropies with c 2 is > 0.06 at the 95% confidence level. When we 
combine CMB anisotropy data with measurements of galaxy clustering, SN-Ia Hubble diagram, and 
other cosmological information, the detection increases to cj is > 0.16 at the same 95% confidence 
level. Future data from Planck, combined with a weak lensing survey such as the one expected 
with DUNE from space, will be able to measure the CNB anisotropy parameter at about 10% accu- 
racy. We discuss the degeneracy between neutrino background ansiotropies and other cosmological 
parameters such as the number of effective neutrinos species and the dark energy equation of state. 



I. INTRODUCTION 

The recent results on the cosmic microwave back- 
ground (CMB) anisotropies from the the five-year 
data of the WMAP mission have once again con- 
firmed the basic predictions of the standard cosmo- 
logical model [l|. In addition to an improved deter- 
mination of several key cosmological parameters such 
as the scalar spectral index and the optical depth to 
reionization with a better control of systematics, the 
new WMAP 5-year data have provided, probably for 
the first time, a clear indication for the presence of an- 
siotropies in the cosmic neutrino background (CNB) 
at more than 95% confidence level. 

The CNB is a clear prediction of the standard cos- 
mological model, though a direct detection of the cos- 
mic neutrinos remain extremely challenging. These 
neutrinos decouples from the primeval plasma at tem- 
perature T ~ 1 MeV before electron-positron annihi- 
lation. A relic neutrino background is expected to- 
day at a temperature of T v — (4/H) 1 / 3 ^ (where 
T 7 = 2.72SK is the CMB blackbody temperature) 
with an energy density of p v ~ 0.58N c gT^. In the 
standard cosmological model N e s = 3.045 and any 
hint for N e g > can be therefore considered as an 
indication for the CNB. 

Recent combined analyses of CMB data with other 
cosmological observables have already provided strik- 
ing evidence for the presence of this neutrino back- 
ground. A combination of WMAP three-year data 
with measurements of Hubble expansion rate, for ex- 
ample, provided the constraint that N c g = 3.7±1.1 at 
95% confidence level [2|, with other analyses finding 
consistent results (see e.g. 3]). The recent WMAP 
5-year data analysis has reported an estimate of N e g, 
but only inluding CMB data in the analysis [lj . Their 
results is then the most conservative, though it is now 
clear that neutrinos arc certainly an important energy 
component of our universe. 

In this paper we reanalyze, in light of the new 
WMAP 5-year data release and their results, some 



further properties of the neutrino background, namely 
the presence of anisotropies in the background and 
not just the effective number of neutrino species. Al- 
though inflationary anisotropies in the CNB at the 
level of ~ 10~ 5 are expected in the standard sce- 
nario, a direct detection of this anisotropy is signif- 
icantly more challenging than a simple detection of 
the relic CNB today. The expected anisotropies in 
the CNB, however, affect the CMB anisotropy angular 
power spectrum at level of ~ 20% through the gravita- 
tional feedback of the neutrino free-streaming damp- 
ing and anisotropic stress contributions [J]. This al- 
lows an indirect detection of the neutrino background 
ansiotropy through its signature in the CMB. 

A way to parameterize the anisotropies in the CNB 
has been introduced in [5J with the "viscosity param- 
eter" c 2 is , which controls the relationship between ve- 
locity/metric shear and anisotropic stress in the CNB. 
In the standard scenario c 2 is = 1/3, anisotropies are 
present in the CNB and approximate the radiative vis- 
cosity of neutrinos. The case c 2 is = 0, on the contrary, 
cuts the Boltzmann hierarchy of CNB perturbations 
at the quadrupole, forcing a perfect fluid solution with 
no CNB anisotropies but only density and velocity 
(pressure) perturbations. Observationally determin- 
ing c 2 is > would therefore provide a strong indica- 
tion for the existence of CNB anisotropies, as argued 
in [6| . This possibility is studied in several papers that 
made use of early CMB data [7| with varying levels of 
indirect detection of the CNB anisotropy. 

In this paper we reanalyze the constraints on neu- 
trino anisotropies in light of the new WMAP 5-year 
data. We combine CMB data with existing large-scale 
structure cosmological information to get an over- 
all estimate on the viscosity parameter and to con- 
straint it better, for the first time, at a confidence 
level greater than 95%. Moreover, we also study how 
future experiments as the Planck satellite CMB mis- 
sion, possibly combined by future weak lensing data, 
will be able to constrain this parameter. Our paper 
is organized as follows: in the next Section, we dis- 



cuss the CMB and large-scale structure data analysis 
and presenr our results based on existing data in Sec- 
tion III. We also forecast the expected errors in Sec- 
tion IV and conclude with a summary of our results 
in Section V. 



II. ANALYSIS 

The method we adopt for this analysis is based 
on the publicly available Markov Chain Monte Carlo 
package cosmomc [Sfl with a convergence diagnostics 
based on the Gelman and Rubin statistic. We sam- 
ple the following eight-dimensional set of cosmological 
parameters, adopting flat priors on them: the phys- 
ical baryon and Cold Dark Matter, 0Jb = Qbh 2 and 
u> c = rich 2 ., the ratio of the sound horizon to the an- 
gular diameter distance at decoupling, S , the scalar 
spectral index ns, the overall normalization of the 
spectrum A at k — 0.05 Mpc -1 , the optical depth to 
reionization, r, the number of massless neutrinos N e g 
and, finally, the viscosity parameter c v i s introduced 
above with the prior c 2 ia < 1/3. Simultaneously we 
also use a cosmic age top-hat prior as 10 Gyr< to < 
20 Gyr. 

Furthermore, we consider purely adiabatic initial 
conditions, we impose flatness and we treat the dark 
energy component as a cosmological constant. We 
include the five-year WMAP data [l| (temperature 
and polarization) with the routine for computing the 
likelihood supplied by the WMAP team. Together 
with the WMAP data we also consider the small- 
scale CMB measurements of ACBAR H, CBI [hJ, 



and BOOMERANG-2K [ll|. In addition to the CMB 



data, we include the Supernovae Legacy Survey data 
[13 |. the real-space power spectrum of galaxies from 
the Sloan galaxy redshift survey (SDSS) [H, and 2dF 
galaxy power spectrum [14J , and, in a separate anal- 
ysis, constraints from the real-space power spectrum 
of red galaxies from the Sloan galaxy redshift survey 
(SDSSlrg) 0. 

Our analysis considers 2 different scenarios: firstly, 
we test for anisotropics in a relativistic background 
comprising N e g = 3 massless neutrino species. Sec- 
ondly, we allow the number of massless neutrinos to 
vary. In both cases, we allow for variations in the 
viscosity parameter c v i s - 

We consider here 3 separate datasets: the first with 
WMAP 5-year only (WMAP5), the second with all 
existing CMB data (ALL) combined with data from 
Sloan red galaxy power spectrum (ALL+SDSSlrg), 
and the third with all CMB information combined 
with the SDSS power spectrum (ALL+SDSS). In the 
last two cases we have also imposed the Guassian 
Hubble Space Telescope prior [l6|: h = 0.72 ± 0.08 
and also a weak big-bang nucleosynthesis prior [17j: 
uj b h 2 = 0.022 ± 0.002 (la). In the computation of 
the CMB spectra we have allowed for the lensing mod- 
ification [lg|. As above for the other parameters we 
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FIG. 1: Joint 2-dimensional posterior probability contour 
plots in the c^ is — n s plane, showing the 68% and 95% con- 
fidence level contours from the WMAP 5-year data alone 
(empty contours), ALL+SDSSlrg dataset (blue contours) 
and ALL+SDSS dataset (red contour). Including the full 
Sloan dataset is necessary to improve the limits placed on 



have adopted flat priors. 



III. RESULTS 

A. Standard background of A c ff = 3 massless 
neutrinos. 

In Figure [T] we plot constraints obtained from our 
analysis in the c 2 is — n s plane in the case of N e g = 3 
massless neutrinos. The WMAP 5-year data alone 
is able to put only a weak constrain on c v i s , due to 
the degeneracy with the spectral index highlighted in 
[lgj l . c 2 - m > 0.06 at 95% confidence level. Including the 
full cosmological dataset and the red galaxy spectrum 
(ALL+SDSSlrg) the degeneracy is partially broken 
enough to obtain cj is > 0.11 (95% confidence level). 
Including the full Sloan + 2dF dataset (ALL+SDSS) 
we obtain a stronger constrain with c 2 is > 0.16 at 
the same 95% confidence level. The difference in 
log-likelihood between the best fit for the standard 
case and the case with no anisotropies (c 2 is — 0) is 
Ax 2 = 15.2. Therefore we can conclude that the case 
where the NB does not have anisotropies above the 
first moment is quite clearly disfavored. 



B. General background of massless neutrinos. 

We now relax the assumption of a standard back- 
ground of relativistic neutrinos by treating N e g as 
a free parameter and thus check the ability of the 
data to put simultaneous constraints both on the 
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scalar spectral index n s together with other cosmolog- 
ical parameters. Several future experiments will focus 
on the observation of weak gravitational lensing with 
large-area galaxy surveys. While weak lensing surveys 
can improve our understanding of neutrino physics 
[2(j |. weak lensing observations help in this analysis 
by providing a way to break certain degeneracies by 
improving the determination of parameters such as 
n s . CMB experiments, however, are more sensitive to 
c^ is , though in CMB alone, this parameter is degen- 
erate with other quantities. It is therefore interesting 
to study how the combination of a planned CMB ex- 
periment, mainly Planck, with weak lensing surveys 
can improve the constraints on the CNB anisotropy 
parameters. 



FIG. 2: Joint 2-dimensional posterior probability contour 
plots in the c^ is — N c g plane, showing the 68% and 95% 
contours from the WMAP5 data alone (empty contours), 
ALL+SDSSlrg dataset (blue contours) and ALL+SDSS 
dataset (red contour). 



A. Fisher matrix 

To explore this we perform a Fisher matrix analysis 
21] where the Fisher matrix is 



background number of neutrinos and on the presence 
of anisotropies in it. In Figure [2] we show the 2- 
dimensional constraints in the plane c 2 is — N e g. 



Dataset 


c 2 . 


A^eff = 


= 3 


N eB 


c?. 


WMAP5 

ALL+SDSSlrg 

ALL+SDSS 


> 0.06 

> 0.11 

> 0.16 


6 4+2-8 

on+1.7 
°- z -i.5 

4 4+2-1 


> 0.07 

> 0.11 

> 0.15 



TABLE I: The 95% confidence level limits on N cS and 
lower limit at 95% c.l. on c 2 is for WMAP5 data alone, 
ALL+SDSSlrg and ALL+SDSS (see text for details). 



As we can see, the WMAP data alone can pro- 
vide a weak indication for for a background of rel- 
ativistic neutrinos, giving 3.2 < N v < 9.2 at 95% 
c.l., and it is able to put also a weak constraint on 
c v i s , giving c 2 ia > 0.07 at 95% c.l.. Considering the 
ALL+SDSSlrg dataset we obtain 1.7 < N v < 4.9 at 
95% c.l. and c^, > 0.11 at 95% c.l.. Including the full 



Sloan dataset increases the constraining on c 
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c z Vls > 0.15 
2.6 < N v < 6.5 



(95% c.l., 1-tail) 
(95% c.L, 2-tails). 



(1) 
(2) 



The difference in log-likelihood between the best fit 
for the standard case and the case with no anisotropies 
(4 S = 0) is A X 2 = 19.4. 



IV. FORECAST 

In this section we investigate how future experiment 
could improve the constraints on the viscosity parame- 
ter cj is , the number of relativistic species N a g and the 
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where L is the likelihood function of the set of param- 
eters pi . Here we adopt a cosmological model with ten 
parameters. We combine the Fisher matrix of Planck 
with the Fisher matrix of three different weak lensing 
surveys to study the dependence of the constraints on 
the survey parameters since weak lensing surveys are 
still under development. According to the Cramer- 
Rao inequality for unbiased estimators this is the best 
statistical error that one can obtain for the generic 
parameter p a : 



by 



<J Pa > ^{F- X ) aa . 
For a CMB experiment the Fisher matrix is g 
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where the C" are the well known power spectra 
for the temperature (TT), temperature-polarization 
(TE), E-mode polarization (EE), and B-mode polar- 
ization (i and j run over TT, EE,TE, BB) and Cov; 
is the spectra covariance matrix. We use the exper- 
imental configuration of the Planck satellite to com- 
pute the Fisher matrix for CMB anisotropies. The 
specifications for this experiment are listed in Table 
ILT1 We use / max = 1500 to calculate the sum in ([3]). 

The future weak lensing surveys allow the possibil- 
ity to "tomographycally" reconstruct the distribution 
of foreground matter responsible for distortions in the 
images of distant background galaxies [22|]. One can 
therefore divide the survey into various redshift bins 
and reconstruct the convergence weak lensing power 
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We choose the probability p(z p h\z) to be Gaussian: 



TABLE II: Experimental specifications for the Planck 
satellite. Channel frequency is given in GHz, FWHM in 
arcminutes, and noise in 10 -6 . 



spectrum Py (£) at multipole £\ the subscripts i and 
j are referred to the redshift bins. The convergence 
power spectrum contains the non linear matter power 
spectrum P n i at redshift z, obtained correcting the 
linear power spectrum P(fc, z). 

With these assumptions the Fisher matrix for weak 
lensing is given by [23 ]: 
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where A£ is the step used for £ and: 



Cjk = P]k + 5 J k{ll nt )n j 1 1 
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(7) 



is the lensing covariance. Here, we ignore non- 
Gaussian contribution to the covariance [? ] as we 
are interested in obtaining an estimate on the ability 
of these experiments to measure CNB parameters. 

In the last expression 7i„ t is the rms intrinsic shear 
(and we assume (jfnt) 1 ^ 2 — 0.22 consistent with high 
fidelity space-based imaging data and future ground- 
based data) and rij is the number of galaxies per stera- 
dian belonging to jth bin: 



3600d 
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fj, 



(8) 



where d is the number of galaxies per square ar- 
cminute and fj is the fraction of sources belonging 
to the jth bin. In our analysis we use £ in the range 
10 < £ < 2000 , so that we are well within the lin- 
ear to mildly non-linear regime of fluctuations and 
not to bias our results with uncertainties in predic- 
tions related to non- linear regime at £ > 2000.. We 
considered three different future weak lensing surveys, 
PanSTARRS [H, DES [H and DUNE [27[ experi- 
ments, with different sky coverage (fsky), median red- 
shift zq and number of sources per square arcminute d. 
We assume a radial distribution function of galaxies 
given by D(z) = z 2 exp[— (z/zo) 1 ' 5 ] where zq depends 
on the survey considered. The experimental charac- 
teristics of these survey are listed in table IIII1 

Here, we account for uncertainties in the measured 
photometric redshifts of the galaxies, according to the 
approach of [28|]: if p(z p f l \z) is the probability that a 
galaxy with (real) redshift z is observed at photomet- 
ric redshift z p h then the distribution of galaxies in the 
ith bin is given by Di(z) = J z ^dz ph D(z)p(z ph \z). 



p{z P h\z) 



/2-Ktj z 



exp 



(Zph—Z) 



with a z depending 



on the survey. 

We then combine the Fisher matrix of Planck with 
the Fisher matrix of different weak lensing surveys to 
show how these surveys can improve the constraints 
coming from CMB. Since CMB anisotropies and weak 
lensing of distant galaxies have origin in two distant 
epochs of the history of the Universe we can consider 
them to be independent. Then, the total Fisher ma- 
trix is simply the sum of the Fisher matrix of CMB 
and that of weak lensing: 



v TOT tpCMB 
*<x0 - *a0 



jpWL 



(9) 



experiment 


J sky 


20 


d a z 


DES 


0.13 


0.8 


10 0.05(1 + z) 


PanSTARRS 


0.75 


0.75 


5 0.06(1 + 2) 


DUNE 


0.5 


0.9 


35 0.03(1 + 2) 



TABLE III: Experimental specifications for the three weak 
lensing surveys of our analysis. 



B. Results 

To compute the convergence power spectrum of 
weak lensing we use the code CAMB with the op- 
tion HALOFIT [23| to obtain the non linear matter 
power spectrum under the halo model [301 ]. Then we 
calculate P^ and its derivatives with respect to cos- 
mological parameters to construct the Fisher matrix. 
We consider the redshift range < z < 3 and divide 
it in 5 redshift bins of equal size. 

For the CMB Fisher matrix we use a modified ver- 
sion of CAMB, to allow variation in c^ is . The deriva- 
tives that appear in §5§ and © are calculated from 
the target model. We use a set of nine cosmological 
parameters whose target values are: ili,h 2 — 0.0223, 
n c h 2 = 0.114, fl A = 0.7, n s = 1.0, r = 0.084, w = -1, 
A s = 2.4 • 10~ 9 , N cS = 3.04 and c^ is = 1/3. We as- 
sume a flat universe imposing the flatness condition: 



h = V(n c h 2 + n b h 2 )/(i - n A ), 



(10) 



In table [IV] we report the la uncertainties for the 
parameters of our model from Fisher analysis. These 
results shows that Planck alone will achieve better 
constraints on the cosmological parameters compared 
to weak lensing surveys, except for the dark energy 
parameters fl\ and w which are very sensitive to the 
tomograhic information coming from weak lensing. 
When combining CMB experiments with weak lens- 
ing surveys, the improvement in la uncertainties is 
of more than one order of magnitude on Qa and w 
and of a factor 1.5-2 on the others parameters. This 



holds even for the viscosity parameter c^ is (although 
this parameter has no effect on the matter power spec- 
trum and hence on the convergence power spectrum 
of weak lensing) because of the degeneracies between 
the parameters of the model, such as the degeneracies 

n s- c vis and TVcff-Cvis- 

The combination of CMB experiments with weak 
lensing surveys can put strong constraints on w, 
though CMB alone can contrain w only weakly. More- 
over Figure [5] shows that there is not any relevant de- 
generacy between w and cj is and a measure of dark 
energy equation of state implies nothing for possible 
anisotropics in the neutrino background. 
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TABLE IV: la errors on the cosmological parameters of 
our model from Planck and weak lensing surveys alone 
(top) and from Planck combined with each survey (bot- 
tom). 



In Figure [31 Figure 0] and Figure [5] we show the joint 
constraints on c^ is with n s , N e g and w, from Planck 
alone and from Planck combined with the three weak 
lensing survey. 

To obtain the confidence contours for any couple of 
parameters, we take the elements of the inverse Fisher 
matrix {F~ 1 ) a p that corresponds to those parameters. 
This gives the correlation matrix of that couple of pa- 
rameters and the eigenvalues and eigenvectors of this 
correlation matrix give the orientation and the size of 
the ellipsoid of confidence. This procedure is equiva- 
lent to marginalize on the remaining parameters (see 
also H3). 



In performing the Fisher analysis we are assuming 
that the likelihood is Gaussian and that the maxi- 
mum of the likelihood is in the target model. The 
likelihood for the parameter c^ is could not be exactly 
Gaussian, at least for c^ is > 1/3. However expand- 
ing the likelihood in Taylor series around the target 
model we can see that the likelihood function is ap- 
proximately Gaussian for c^ is < 1/3, (i.e. it is domi- 
nated from quadratic terms), even if it could be not 
Gaussian globally. In fact usually L drops rapidly so 
that we can neglect the non-Gaussianity of the like- 
lihhood [32| and the constraints from Fisher analysis 
can be considered a good approximation. 

The most considerable result of this analysis is 
that it is possible to detect neutrino background 
anisotropies even from Planck alone with 0.33 > 
c^ is > 0. The combination of Planck experiment with 
weak lensing surveys can decrease the statistical un- 
certainty on c^ is about of a 35% (for Planck+Dune), 
while we find the combinations Planck+DES and 
Planck+PanSTARRS give similar constraints. 
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FIG. 3: Constraints for c 2 iB and n s from Planck alone 
and from Planck combined with DUNE experiment (top) 
and from Planck combined with others redshift surveys 
(bottom). 
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FIG. 4: Same as Figure [3] but for c^, ia and N e g 



result. 

In this respect we have performed a forecast for fu- 
ture CMB and weak lensing missions. We have found 
that a combination of DUNE and Planck will be able 
to yield a measurement of the c v - ls parameter at ~ 10% 
level. Moreover, we studied the possible degeneracies 
with other cosmological parameters. Including uncer- 
tainties on c v i s will double the error bars on N e g while 
the constraints on the equation of state w will remain 
pratically unaffected. While the standard model pre- 
dicts c^ is = 1/3 several cosmological scenarios, from 
interacting neutrinos to early dark energy, can be con- 
sidered that could bring a deviation in the measured 
value. A better detection of anisotropies in the neu- 
trino background could therefore provide an useful 
test for those models and possibly indicate the pres- 
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V. CONCLUSIONS 



FIG. 5: Constraints on c^ is and w from Planck combined 
with the three weak lensing surveys. 



In this paper we have re-analyzed the status of neu- 
trino anisotropies in light of the new WMAP five-year 
data. We found that while the WMAP 5-year data 
alone are unable to provide significant evidence for 
those anisotropies, combination of CMB data with 
current large-scale structure cosmological data yields 
a detection at more than 2cr confidence level. Future 
cosmological data are certaily needed to confirm this 



ence of new physics. 
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